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Abstract. Totally asymmetric simple exclusion processes, consisting of two coupled 
parallel lattice chains with particles interacting with hard-core exclusion and moving 
along the channels and between them, are considered. In the limit of strong coupling 
between the channels, the particle currents, density profiles and a phase diagram 
are calculated exactly by mapping the system into an effective one-channel totally 
asymmetric exclusion model. For intermediate couplings, a simple approximate 
theory, that describes the particle dynamics in vertical clusters of two corresponding 
parallel sites exactly and neglects the correlations between different vertical clusters, 
is developed. It is found that, similarly to the case of one-channel totally asymmetric 
simple exclusion processes, there are three stationary state phases, although the phase 
boundaries and stationary properties strongly depend on inter-channel coupling. An 
extensive computer Monte Carlo simulations fully support the theoretical predictions. 
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1. Introduction 

In recent years asymmetric simple exclusion processes (ASEPs) have become an 
important tool of investigation for many processes in chemistry, physics and biology 
[H Ej. ASEPs have been applied successfully to understand the kinetics of 
biopolymerization [3 , polymer dynamics in dense medium [3], diffusion through 
membrane channels jS] , gel electrophoresis [H] , dynamics of motor proteins moving along 
rigid filaments [Zj , and the kinetics of synthesis of proteins |H1 EH HH| . 

ASEPs are one-dimensional lattice models where particles interact only with hard- 
core exclusion potential. Each lattice site can be occupied by a particle or it can be 
empty. In the simplest totally asymmetric simple exclusion process (TASEP), where 
the particles can only move in one direction, the dynamic rules are the following. The 
particle at site % can jump forward to the site i + 1 with the rate 1 if the target site 
is empty. The particle can enter the lattice with the rate a, provided the first site is 
unoccupied, and it can also leave the system with the rate (3. Although these dynamic 
rules are very simple, they lead to a very rich and complex dynamic phase behavior. 
There are non-equilibrium phase transitions between the stationary states of the system, 
induced by boundary processes, that have no analogs in equilibrium systems [T| l2*| lllj. 

The coupling of ASEPs with different equilibrium and non-equilibrium processes 
have led to many unusual and unexpected phenomena. An introduction of a single 
irreversible detachment process to the bulk of the system have resulted in significant 
changes in phase diagram The coupling of ASEPs with equilibrium Langmuir 

kinetics at each site produced unusual phenomena of localized density shocks [121 El 
EH E3 ■ However, the complexity of phase transitions in ASEPs can be reasonable well 
explained and understood by applying a phenomenological domain- wall theory |2*| lllj. 

Most investigations of ASEPs concentrate on one-channel systems where particles 
can only move along one lattice chain. However, the description of many real phenomena 
would be more adequate if parallel-chain asymmetric exclusion processes are used. 
For example, motor proteins kinesins can move along the parallel protofilaments 
of microtubules, and there are no restrictions for them to jump between these 
protofilaments 17 . Parallel-chains ASEPs have been considered earlier ^SJ El EH]- 
However, in these investigations the coupling between different chains was indirect, 
i.e., hopping between the chains was forbidden. The aim of the present paper is to 
investigate two-chain asymmetric exclusion processes, where particles can move between 
the lattice channels. We investigate this system by using a simple approximate model 
for intermediate couplings and by mapping it to the exactly solved one-channel TASEP 
in the limits of strong coupling and no coupling. In addition, extensive computer Monte 
Carlo simulations are performed in order to validate theoretical predictions. 

The paper is organized as follows. In section 2 we give a detailed description of the 
model, discuss known results for one-channel ASEPs, solve exactly the stationary state 
properties of the system in the limit of strong coupling, and develop an approximate 
theory for two-channel asymmetric exclusion processes for intermediate couplings. Then 
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Figure 1. Schematic view of the model for two-channel TASEP. Allowed transitions 
are shown by arrows. Inter-channel hopping rates are equal to w. The transition rates 
to move along the channel are 1 — w, if there is no particle at the same site in another 
channel, otherwise it is equal to 1. Entrance and exit rates are a and /3 (or /3(1 — w) 
if the last cluster is half-empty), respectively, and are the same in both channels. 



in section 3 we present and discuss Monte Carlo simulations results. Finally, we 
summarize and conclude in section 4. 



2. Theoretical Description of Two-Channel TASEP 

2.1. Model 

We consider identical particles moving on two parallel one-dimensional lattices, each 
lattice has L sites, as shown in Fig. ^ Every site on both lattices can be either empty 
or it can be occupied by no more than one particle. In the bulk of the system the 
dynamic rules are the following. A particle at site i can hop up or down to the same 
site i on the other lattice chain with the rate < w < 1, if that site is empty. The 
particle can also move from left to right along the same lattice to site i + if this site is 
available. However, this transition rate depends on occupation of site i at neighboring 
channel. If there is no particle at that site, the rate is equal to 1 — w, otherwise the 
particle jumps with the rate 1 (see Fig. Q). It means that the full transition rate of 
leaving site i, to go forward or up/down, is always equal to 1. Particles can also enter 
the system with the rate a if any of the first sites in either channel are empty. When a 
particle reaches site L, it can exit with the rate /3, when both last sites are occupied, or 
with the rate (3(1 — w) if there is no vertical neighbor at the other lattice chain. 
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When a vertical transition rate w = 0, we have two uncoupled single-channel totally 
asymmetric simple exclusion processes for which an exact description of phase diagram, 
density profiles and particle currents is known PJ E] . In this case there are three steady- 
state phases, determined by the processes at the boundaries or in the bulk of the system. 
For a > 1/2 and (3 > 1/2, the particle dynamics is specified by the processes in the 
bulk of the system, and we have a maximal- current phase with the following stationary 
current and bulk density, 

JhlC = 1/4, Pbulk,MC = 1/2. (1) 

The entrance into the system determines the overall particle dynamics for a < (3 and 
a < 1/2. In this case the system can be found in a low-density phase with the current 
and bulk density given by 

Jld = «(1 - a), p bu ik,LD = ct. (2) 

When the exit from the system limits the overall dynamics, which takes place for j3 < a 
and (3 < 1/2, the system exists in a high-density phase. In this case the current and 
bulk densities are equal to 

Jhd = (3(1 - (3), Pbuik,HD = 1 - /3. (3) 

For vertical transition rates w > 0, the particle dynamics in both channels depend 
on each other, and any successful theoretical description should be able to account for 
this coupling. Let us consider a cluster of two vertical sites i. There are four possible 
states for this cluster, as shown in Fig. EJP^ is defined as a probability to find a cluster 
with both lattice sites occupied. Then P$ and P ^ are the probabilities that only lower 
or upper channel site is occupied (see Fig. |2J). Finally, P q is a probability that both 
sites in the cluster are empty. The normalization condition for these probabilities gives 
us 

P$?+P$!+P$!+P!8 = l. (4) 

In the bulk of the system at stationary state, it is reasonable to expect that these 
probabilities are independent of the position of the vertical cluster due to translational 
symmetry. In what follows, we omit the superscript % for bulk values of these 
probabilities since only steady-state processes will be considered. Also, because both 
channels are equivalent, we expect that p/q = Pqi- Then the bulk density at each 
lattice chain, which is also independent of the position on the lattice in stationary state 
limit, is given by 

p = P u + Pio. (5) 

Thus the overall dynamics of the system can be fully described in terms of these 
probabilities. 
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Figure 2. Four possible states for a vertical cluster of lattice sites. Pn, P%o, Pqi and 
Poo are the corresponding probabilities for each state. 



2.2. Strong- Coupling Limit 

In the limit of strong coupling (w = 1), the dynamics of the system simplifies 
significantly, because at large times any vertical cluster cannot exist in a configuration, 
where both sites are empty, i.e., Poo — 0. This can be seen from the following arguments. 
A vertical cluster in state {00} could only be obtained by moving particle along the 
corresponding lattice chain if the previous state of the same cluster was {10} or {01}. 
However, the rate for this transition in the limit of strong coupling is equal to 1 — w — 0, 
and the configuration {00} can never be reached for any vertical cluster. Thus, in this 
system there are only two types of clusters: fully filled and half-filled. 

One can think of filled clusters as new effective "particles" and half-filled clusters 
then can be viewed as new effective "holes." The two-channel TASEP in this limit can 
be mapped into the one-channel totally asymmetric exclusion process with the effective 
entrance rate a and the effective exit rate 2(3. In the bulk of the system an effective 
new particle jumps to the right with rate 1, if this move is allowed. The factor 2 in the 
effective exit rate is due to the fact that there are two ways of leaving the system from 
the filled cluster at last site, namely, by moving particles from the upper or lower lattice 
channels. 

The exact density profiles, particle currents and phase diagram for this effective 
one-channel TASEP system are known |UI2j, and it allows us to calculate exactly the 
properties of the original two-channel TASEP in the limit of strong coupling. The 
steady-state particle current of the effective one-channel model J* is related to the 
particle current per channel J of the original two-channel model in the following way, 

j = r/2, (6) 

while using Eq. |S] we obtain the corresponding relation for bulk density profiles in 



Two- Channel Exclusion Processes 



6 



two-channel TASEP, 

p = p* + (l-p*)/2 = (l + p*)/2, (7) 

where p* is the effective density for one-channel asymmetric exclusion model. 

Thus we conclude that there are 3 stationary-state phase for two-channel TASEP 
in the limit of strong coupling. When a < 2(3 and a < 1/2, the system can be found in 
the low-density phase with the following properties, 

Pii,ld = a, Pw,ld = (l-a)/2, J LD = a(l-a)/2, p bu ik,LD = (1+a)/ 2 - 

The conditions a > 2(3 and (3 < 1/4 specify the high-density phase. In this case the 
steady-state properties are given by 

P\1,HD — 1 — 2/3, PlO,HD = P, JhD=P(^ — ^P), Pbulk,HD = 1 — /?-(9) 

In the case when bulk dynamics is rate-limiting, we have the maximal- current phase 
with the following parameters, 

P n , MC = l/2, P W , MC = 1/4, Jmc = 1/8, Pwfc.Afc = 3/4. (10) 

<?. Approximate Solutions for Intermediate Couplings 

For intermediate couplings, < w < 1, any vertical cluster can exist in all four 
possible states and we cannot map the two-channel TASEP into the effective one-channel 
model. Some reasonable approximations are needed in order to calculate the steady- 
state properties of the system. 

Assume that a state of a given vertical cluster is independent of states of its 
neighbors. It means that the cluster dynamics is considered in a mean-field description. 
Then the probability of finding the cluster at the position i with both sites occupied 
changes in time as given by 
dP 

= P n (P 10 + P i) + 2(1 - w)P w P 01 - 2P n P 00 - P U {P 10 + P i), (11) 

where the first term corresponds to the probability density flux from sites i — 1, the 
second term describes the changes at both sites i that lead to the formation of the 
cluster {11}, and two negative terms represent the density flux leaving from sites i. 
Applying the symmetry relation P w = P m at large times simplifies this expression into 

= (l-«j)if o -P u Poo. (12) 
This equation, with the help of the normalization condition (0J), yields 



-P n + ^P? 1 + (1- W )p u (l-P u ) 

Pio - : • U-iJ 

1 — w 

Such rather complex relation can take a simpler form for two limiting cases - for two- 
channel TASEP without coupling and in the strong-coupling limit. For w = we obtain 

Pio = JPu-Pii, (14) 
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while for w = 1 it transforms into 

fl. = (15) 

in agreement with exact calculations presented above [see Eqs. (8-10)]. 
The particle current per channel can be written as 

J = [P n + (1 - w)P 10 ] (1 - Pio - Pn). (16) 

In this expression the first multiplier gives the probability to find a vertical cluster in 
configurations, from which the particle can hop to the right, while the second multiplier 
gives the probability that the site ahead is available. Using the relation (|13j) we can 
express the particle current only in terms of one variable - the density of vertical clusters 
Pn, 



■l-w + wP u - J Pull -w + wP ] 



J = JP 11 (1 -w + wP u ) f (17) 

v 1 — w 

Again, simpler relations can be obtained for two limiting cases of zero coupling and 

strong coupling. The particle current in the case of w = is given by 



J = yJPu (1 - yJP n j , (18) 

while for w = 1 we conclude that 

J = Pu(l-P u )/2. (19) 

These relations also agree with exact results discussed above. 

We expect that, similarly to the case of zero inter-channel coupling or in the strong- 
coupling limit, there are three possible stationary phases. The conditions for a maximal- 
current phase can be specified by computing a maximum of the particle current as a 

dJ 
dP n 



function of density of {11} clusters, i.e., -M- = 0. It leads to the following equation, 



l-ur /'.,(! -w + wP n ) + (1 -" 2 ^g l)1,5 -1 + w- 2wP u 



0. (20) 
1 — w 

For any value of inter-channel coupling w this equation can always be solved exactly 
numerically, and we can obtain the value of the density of filled vertical clusters in the 
maximal- current phase. Because this equation can have more than one real solution, we 
should always choose the physically reasonable solution which gives 1/4 < P\\,mc < 1/2. 
Utilizing Eqs. (0), (jlS} and ljT7j) . the stationary properties of this phase can be easily 
determined. For example, for w = 0.5 the calculations show that 

Pll,max - 0-317, Pl0,max ^ 0.280, p bu lk,max ^ 0.597, J max ~ 0.184.(21) 

In the low-density phase the entrance rate a limits the overall particle dynamics, 
and the expression for the current per channel is given by 

J LD = a(l-P u -P 10 ), (22) 

or, using Eq. (JT3J), it can be written as 



1-w + wPn- JPn(l-w + wP n ) 

Jld = a f . (23) 

1 — w 
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Comparing this equation with the expression for the particle current in the bulk of the 
system [see Eq.flTTj)]. it can be shown that 



a =y/P 11 (l-w + wP n ), (24) 

from which the value of density of {11} clusters in the low-density phase can be obtained 
as follows, 

-l + w+ J(l - w) 2 + 4wa 2 

= —„ • < 25) 

Applying again Eq. (JT3J), we derive the expression for the density of {10} clusters, 



10,LD 



a i-V 1 + w lr^ 



+ V - ^ • (26) 

1 — w 2w 

These equations allow us to calculate the bulk density as p = P\\ + Piq, and also, 
substituting Eq. (J2*5|) into Eq. (|23j). the particle current per channel can be expressed 
as 



j _ a 



2a / (a. 
1 - + i/l + 4w 



(27) 

1-w V \l-wj 

When exit dynamics determines the overall behavior of the system we have a high- 
density phase. The particle current per channel is given by 

Jhd = P[Pii + (1-w)P 10 ], (28) 

and, after applying Eq. ([T3]). it transforms into 



J HD = PyJp n {l-w + wPu). (29) 

Again comparing this expression with the particle current in the bulk, as given by 
Eq.([T7|). we derive the density of filled vertical clusters, 

1 - Jl - 4w[3{l - J) 

Pu,hd = 1-P V ^— -. (30) 

2w 

Substitution of it into Eq. (fT3)) allows us to compute the density of {10} clusters, 



P _ Wl- 4^(l-fl 
2w 

Then the bulk density in each channel is equal to phd = P\i,hd + Pio,hd = 1-/3. 
It is interesting to note that the bulk density in this phase is always independent of 
inter-channel coupling and depends only on exit rate, although the densities of filled 
and half-filled clusters strongly depend on w. The expression for the particle current 
can be derived from Eqs. (}2T?j) and (J30J), 



Jhd-~ 



20+y/l- 4wP(l - (5) 



(32) 



There are two types of phase transitions in the system. The boundary between the 
low-density and high-density phase specifies the first-order phase transition with jumps 
in bulk densities, while the transition between the maximal-current and low-density 
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Figure 3. Phase diagrams for two-channel TASEPs for different inter-channel 
couplings. Symbols correspond to Monte Carlo computer simulations, lines are our 
theoretical predictions. 



or high-density phases is continuous. The low-density and high-density phases coexist 
when the particle currents in each phase are equal, i.e., Jhd = Jld- From Eqs. (J27|) 
and (J32j) we obtain, 



a 
2 



2a 



+ \ l+4w 



w 



a 



1 — w 



2 



1-2(3 + yjl - 4wP(l - (3) 



,(33) 



which determines the phase boundary curve in {a,/3} coordinates. One can observe 
that, in contrast to the systems with zero coupling or in the limit of strong coupling, 
the phase boundary here is not a linear function of parameters. Similar arguments can 
also be used to find the boundaries between the low-density or high-density and the 
maximal- current phases. Theoretically calculated phase diagrams for different inter- 
channel couplings are presented in Fig. H3 



3. Monte-Carlo Simulations and Discussions 



Our approximate theory treats exactly the particle dynamics inside vertical clusters, 
however, interactions between the clusters are accounted for in an approximate mean- 
field fashion. Theoretical predictions from this approximate method agree with exact 
results in the limiting cases of zero inter-channel coupling and strong coupling. To check 
the overall validity of our approach for intermediate couplings we performed extensive 
Monte Carlo computer simulations. 

Our theoretical arguments are valid only in the thermodynamic limit, i.e., L — > oo. 
However, in our simulations the lattice of size L = 100 was used, and we checked that 
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Figure 4. Density profiles for inter-channel coupling w = 0.5: a) in the high-density 
phase with a = 0.8 and f3 — 0.2; b) in the low-density phase with a = 0.2 and f3 = 0.8; 
c) in the maximal-current phase with a = 0.8 and (5 = 0.8; and d) at the boundary 
between the low-density and the high-density phases for a = 0.2 and j3 = 0.2. Symbols 
correspond to Monte Carlo computer simulations results, lines describe theoretically 
calculated bulk values for Piq, Pn and density per channel, respectively. Error bars, 
determined from standard deviations for simulations, are smaller than the size of the 
symbols. 



for larger lattices there is no difference with reported results. The density profiles and 
particle currents in our simulations were calculated by averaging over 10 8 — 10 10 Monte 
Carlo steps, although first 5% of total number of steps were neglected to ensure that the 
system has reached the stationary state. Boundaries in phase diagrams were determined 
by comparing density profiles and changes in the particle current. 

Phase diagrams and density profiles calculated from Monte Carlo simulations are 
presented in Figs. El and |U The agreement with theoretically calculated values is 
excellent. It should be noted that the phase boundary between the low-density and 
the high-density phases for intermediate coupling (for < w < 1) deviates slightly from 
linear relation, in agreement with theoretical predictions. That can be seen by analyzing 
the slope of the curve. It equals to 1 at small a and (3 and slowly decreases for higher 
values of parameters, as shown in Fig. El The overall effect of inter-channel couplings 
is the decrease in a phase volume for the high-density phases. Note, that although we 
presented only theoretical bulk values for densities of clusters and particle densities (see 
Fig. HJ), the approximate theory allows us also to calculate the full density profiles. 

It is interesting also to analyze the effect of inter-channel particle transitions on the 
dynamics of two-channel TASEPs. As shown in Fig. 03 the increase in inter-channel 
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coupling generally decreases the particle current per channel, while the values of bulk 
densities do not change (for the high-density phase) or go up. The strongest effect is 
for the maximal-current and the low-density phases, while the influence of coupling on 
the high-density phase is quite minimal. That can be easily understood by recalling 
what processes determine the different steady-state phases. The entrance processes and 
bulk processes are strongly affected by inter-channel particle transitions, while the exit 
processes do not depend much on this coupling. The inter-channel coupling has also a 
peculiar effect on a triple point in phase diagram, where all three phases coexist - see 
Fig. |U1 The increasing inter-channel current decreases /3-coordinate of the triple point, 
while the effect on a-coordinate is non-monotonous. 

4. Summary and Conclusions 

We investigated two-channel totally asymmetric simple exclusion processes with the 
possibility of particle transitions between the lattice chains. When there is no coupling, 
w — 0, the system can be viewed as two independent one-channel TASEPs, for which 
exact solutions are known. In the limit of strong coupling, w — 1, the exact description 
of particle dynamics is obtained by mapping the two-channel system into an effective 
one-channel TASEP with known stationary properties. 

The two-channel TASEPs with intermediate couplings are analyzed with the 
help of an approximate theoretical approach, in which all stationary properties are 
obtained analytically. In our approach the particle dynamics inside vertical clusters of 
corresponding lattice sites is considered exactly, while the correlations between different 
clusters are neglected. The results of this approximate method are in excellent agreement 
with computer Monte Carlo simulations. Our theoretical and computational results 
indicate that the inter-channel currents have a strong effect on steady-state properties of 
the system. Increasing the coupling lowers the particle current per channel, increases the 
bulk values of particle densities and shifts significantly the position of phase boundaries, 
although the overall topology of the phase diagram is preserved. 

There are several extensions of the original two-channel TASEP that will be 
interesting to investigate. In this paper we considered only the case with symmetric 
coupling, i.e., for any particle the probability of jumping to a different lattice chain is 
independent of the channel. It will be interesting to analyze the asymmetric coupling, 
where particles in the first channel could hop to the second channel with the rate 
Wi, while the opposite motion has the rate of u>2, and w\ ^ u>2- Another more 
interesting, although more complex, case is the problem of non-uniform couplings 
between the channels. It seems reasonable to suggest that our method of combining 
simple approximate theory with computer Monte Carlo simulations is a promising 
approach to study these complex non-equilibrium one- dimensional problems. 



Two- Channel Exclusion Processes 



12 



[a] 



0.05 

°. 
1 




• 


LD sim 




LD theor 


■ 


HD sim 




HD theor 


T 


MC sim 




MC theor 



0.2<h 







[b] 



± 



0.2 



0.4 0.6 

w 



0.8 




• 


LD sim 






■ LD theor 




■ 


HD sim 






- HD theor 




▲ 


MC sim 






- MC theor 





0.2 



0.4 



0.6 



0.8 



w 



Figure 5. Particle currents per channel and bulk densities as a function of inter- 
channel coupling. The low-density phase is specified by a = 0.2 and (3 = 0.8; the 
high-density phase is given by a — 0.8 and j3 — 0.2; and the maximal-current phase is 
for a — 0.8 and (3 = 0.8. Symbols correspond to Monte Carlo computer simulations, 
lines are our theoretical predictions. 
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Figure 6. The coordinates of triple points for different inter-channel couplings. Monte 
Carlo computer simulations results coincide with the theoretical calculations, and thus 
are not shown. 
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